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ABSTRACT 



In October 2007 a hard X-ray burst was detected by the INTEGRAL satellite from a direction consistent with the position of AX 
J1818.8 1559, an X-ray source at low Galactic latitude discovered with the ASCA satellite in 1996-1999. The short duration (0.8 
s) and soft spectrum (power law photon index of 3.0±0.2) of the burst in the 20-100 keV range are typical of Soft Gamma-ray 
Repeaters and Anomalous X-ray Pulsars. We report on the results of an observation of AX J 1 8 1 8.8 — 1559 obtained with the Suzaku 
satellite in October 2011. The source spectrum, a power law with photon index ~1.5, and flux (~ 2 x 10~ 12 erg cm" 2 s , 2-10 
keV) do not show significant variations with respect to the values derived from archival data of various satellites (ROSAT, XMM- 
Newton, Chandra, Swift) obtained from 1993 to 2011. We discuss possible interpretations for AX J1818.8 — 1559 and, based on 
its association with the INTEGRAL burst, we propose it as a new member of the small class of magnetar candidates. 
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1. Introduction 



Although early suggested (e.g. lPaci"nilll967l) . the relevance of 
magnetic energy in powering the emission from neutron stars has 
been recognized only recently, with the discovery of anomalous 
X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs). 
These are spinning-down, isolated neutron stars with luminos- 
ity in the soft and hard X-rays typically larger than their rota- 
tional energy loss and characterized by the emission of power- 
ful bursts and flares. AXPs and SGRs were historically divided 
in two classes, but many observations in dicate that there a re no 
substantial differences betwen them (see Mereghetti 2008 for a 
review ). Even though alterna t ive interpretations have been pro- 
posed dChatterjee et : alj|2000l: |Alpdl200ll ICeall2006t |Xu|[2007; 
lHorvathll2007l: buved et alj|2007t lErtan et alj|2009f the model 
involving magnetars, i.e. highly magnetized neutron stars, is the 
one th at currently better explains the properties of AXPs and 
SGRs dDuncan & Thompsonlll992l:rThompson & Duncanlll995t 
iThompson et al]l2002l) . According to the magnetar model, the 
emission from these sources is ultimately powered by the energy 
stored in their strong magnetic fields, with B ~ 10 13 -10 14 G in 
the magnetosphere, and possibly reaching even larger values in 
the neutron star interior. 

The originally small sample of AXPs/SGRs has increased 
significantly in the latest years, mainly thanks to the discov- 
ery of transient sources, which remained undetected (or un- 
noticed) until they underwent bright outbu rsts, often accompa - 
nied by the emission of short bursts (Re a & Esposito! 12011). 
In a few other cases, new magnetar candidates were found by 
thorough investigation of persistent X-ray sources with peculiar 
prope rties and/or possible associations with supernova rem nants 
(e.g. lGelfand & Gaensled2007tlHalpern & Gotthelfl2010blah. or 
by the discovery of a peculiar radio pulsar dLevin et alj l2010t 
And erson etai1l2012l) . 
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Many of the new additions to the magnetar family have 
shown that these objects present a variety of properties unsus- 
pected from the original members of the class and which connect 
them with other classes of isolated neutron stars. The radio de- 
tectio ns ofXTEJ1810-197and IE 1547-54 dCamilo et alJ 2006. 
120071) and the v ery small spin-down rate of SGR 04 18+5729 
dRea et alJl2010h . implying a dipole field close to that of normal 
pulsars, are just a few examples. 

The X-ray source AX J18 18.8— 1559 was discovered during 
a surv ey of the Galactic pl ane carried out with ASCA in 1996- 
1999 (Sugizak ietaill200lh . but it did not attract particular at- 
tention until INTEGRAL detected a short burst from its direction 
in October 2007. Here we report on a detailed re-analysis of the 
INTEGRAL burst and on the results of a long X-ray observa- 
tion of AX Jl 8 18.8 — 1559 carried out with the Suzaku satellite 
in October 201 1. We also present archival X-ray data, that were 
mostly unpublished, to provide an exhaustive description of the 
X-ray properties of AX J1818.8 — 1559, and discuss its possible 
inclusion in the class of magnetar candidates. 

2. The burst of 2007 October 17 

A short burst at low Galactic latitude (b=-0.26), was discov- 
ered and locali zed in real time with th e INTEGRAL Burst Alert 
System (IBAS; Mereghe tti et alj|2003l) on 2007 October 17 . The 
burst, reported as GRB 071017 dMereghetti et af] 120071). wa s 
found in the data of the ISGRI detector dLebrun et alJ 120031) . 
which is the low-energy imag er of the IBIS code d-mask in- 
strument on board INTEGRAL dUbertini et alJl2003h . The burst 
light curves in different energy ranges (Fig. Q]) show that, un- 
like the typical gam ma-ray bursts (GRBs) observed by IBIS 
dVianello et alJ2 009). this event was not detected above 100 keV. 
The 15^-0 keV light curve is characterized by a single pulse 
lasting ~0.8 s, with a fast rise time starting at 00:58:08.55 UT 
followed by a slower decay. 
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Table 1. Log of the observations. 



Satellite 


Obs.Id 


date 


Duration (ks) 


Notes 


ROSAT 


rp50311nOO 


1993 September 12 


6.7 


PSPC 


XMM-Newton 


0152834501 


2003 March 28 


10.0 


EPIC 


Swift 


00020057001 


2007 October 19 


10.4 


XRT 


Chandra 


7617 


2007 October 21-22 


5.1 


ACIS-S TE 


Chandra 


7616 


2007 October 22 


39.9 


ACIS-S CC 


Chandra 


9609 


2008 February 19 


1.6 


ACIS 


Chandra 


9610 


2008 September 25 


1.2 


HRC-I, 12' off-axis 


Suzaku 


406074010 


2011 October 15-18 


80 


XIS 


Swift 


00044118001 


2012 March 4 


0.5 


XRT 



We performed a spectral analysis by extracting the ISGRI 
data of the first 0.2 s, corresponding to the brightest part of the 
burst. The results confirm that, as suggested by the light curves 
in different energyt bands, the burst was rather soft. Its spectrum 
is equally well fit by a power law with photon index 3.0 ± 0.2 
or by a thermal bremsstrahlung with temperature kT=25±3 keV 
(Fig.|2]i. In both cases the average flux in the 20-100keV energy 
range is ~ 2 x 10~ 7 erg cm -2 s _1 . 

The burst localization that was obtained in real time by the 
automatic IBAS software had an error of 2.9' and was consis- 
tent with the position of AX J1818. 8-1559. We reanalyzed the 
ISGRI data using the most recent versions of the satellite atti- 
tude information and imaging software. Using the data in the 
15-50 keV, where the signal to noise is maximum, we obtained 
a 11.5(7 detection of the burst at the coordinates R.A.=274.731, 
Dec.=- 15.993 (J2000), with an uncertainty of 1.8' (90% confi- 
dence level radius). The Corresponding galactic coordinates are 
1=15.04, b=-0.26. The INTEGRAL error region is indicated by 
the circle in Fig. [3] 

3. XMM-Newt on archival data 

The region of sky in the direction of GRB 071017 was ob- 
served for about 10 ks on 2003 March 28, during a survey of 
the Galactic plane carried out with the XMM-Newton satellite. 
Using the Science Analysis Software (SAS VI 1) we reprocessed 
the data o f the EPIC instr ument (0.2-12 ke V), which consists of 
one pn dTurner et al.ll2001l) and two MOS dStruder et alj|2001l) 
cameras. A few X-ray sources are detected in the EPIC im- 
ages, but only one, at an off-axis angle of 4', lies in the er- 
ror circle of the INTEGRAL burst (Fig. 0). Its coordinates are 
R.A. = 18 h 18 m 51!5, Dec. = -15°59'23" (J2000), with an 
uncertainty of 2". Based on the positional coincidence, as well 
as on the consistent spectral properties (see next paragraph), 
we identify the XMM -Newton source with AX J1818. 8-1559 
(ISugizaki et al.ll200lh . 

No other sources are detected with XMM-Newton inside the 
error circle of the burst, down to a flux limit of a few 10~ 14 erg 
cm -2 s _1 . A few other sources are visible in the EPIC image. 
A faint source is located 1' from the southern border of the error 
circle, at R.A. = 18 h 18 m 57:9, Dec.=-16° 02' 08" (J2000). Its 
flux is about a factor 20 smaller than that of AX J18 18.8- 1559. 
A brighter source, at R.A. = 18 h 19 m 13?6, Dec.=-16° 01' 20" 
(J200 0), can be identified with the ASCA source AX J181915- 
1527 (ISugizaki et al.ll2001l) . 

To study the spectrum of AX J18 18.8— 1559, we extracted 
the pn and MOS source counts from a circle of 40" radius and 
the background spectra from source free regions of the same 
CCD as the source. After rebinning to have at least 30 counts in 
each spectral bin, we fitted simultaneously the pn and two MOS 



spectra. A good fit was obtained with a power law with photon 
index T = 1.2 ± 0.2, absorption /V H = (3.6 ± 0.5) x 10 22 
cm~ 2 , and absorbecQ flux 1.4 x 10~ 12 erg cm~ 2 s _1 in the 2- 
10 keV energy range. These results are consistent with those 
we obtained in a quick look analysis of the sa me data per- 
forme d soon after the INTEGRAL GRB detection (Tiengc Tet al.l 
2007). The source distance is unknown, but considering its lo- 
cation in the direction of the Galactic center region and the 
large absorption, we assume in the following a reference value 
of 10 kpc. The 0.5-10 keV luminosity, corrected for the ab- 
sorption, is ^2.4 x 10 34 erg s _1 . An equally good fit is ob- 
tained with a blackbody of temperature kTs b = 1.9 ± 0.1 keV 
and emission radius Rbb = 0.12 ± 0.05 km. In this case the 
absorption is Nn = (1.6 ± 0.3) x 10 22 cm -2 . For compari- 
son, th e total Galactic absorption in this direction is 1.2xl0 22 
cm~ 2 (Kalberla et al. 2005). All the spectral results are summa- 
rized in Table |2] The values obtained with the power law model 
are very similar to those measured with ASCA in 1996-1999 
(ISugizaki et al.ll2001l) . 

We searched for pulsations in the EPIC data with Fast 
Fourier Transform (FFT) and Rayleigh test methods without 
finding any statistically significant signal. However, this result 
is not particularly constraining due to the relatively small num- 
ber of counts. The 3cr upper limits on the pulsed fraction, A, in 
the 1-10 keV energy band are of 50% for periods in the 0.3-1000 
s range and of ^50-70% between 145 and 300 msB- 

4. Suzaku observation 

We performed a Suzaku observation of AX J 1 8 1 8 . 8 — 1559 from 
2011 October 15 at 13:17 UT to October 18 at 03:32 UT. This 
provided about 80 ks of exposure in the XIS instrument. The XIS 
consists of four telescopes with a spatial resolution of about 2' 
coupled to four CCD cameras o perating in the 0.2-12 keV en- 
ergy range (Kovama et al. 2007). At the time of our observation 
only the two front illuminated XIS0 and XIS3, and the back- 
illuminated XIS1 were operating. The three detectors were used 
in 1/8 window mode, which gives a time resolution of 1 s at the 
expense of a reduced field of view. 

For the spectral analysis we screened the data using standard 
criteria. In particular, we considered only events with GRADE = 
0,2^4,6 and excluded time intervals within AT = 436 s from the 
South Atlantic Anomaly (SAA) and pointed at elevation angles 

We use the PHABS absorpt ion model in XSPEC with abundances 
from lAnders & Grevessd (1989). 

2 Here and in the following we assume a sinusoidal modulation of the 
flux, F(t)=F [l+Asin(Slt)] and define the pulsed fraction as the sinusoid 
amplitude A. All the upper limits on A are computed as described in 
IIsrael&Stell3H996l) 
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Fig. 1. Light curves of the 2007 October 17 burst measured with ISGRI 
in three energy ranges. The bin size is 50 ms. 
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Fig. 2. ISGRI spectrum of the 2007 October 17 burst. Top panel: data 
and best fit power law model. Bottom panel: residuals in units of sigma. 



smaller than 5° or 20° above the night or day Earth rim, respec- 
tively. Source and background spectra were extracted from cir- 
cles of 1.1' radius and rebinned so as to have a minimum of 20 
counts per spectral bin. Slightly less stringent criteria were used 
for the timing analysis (AT = 60 s, 10° day Earth minimum el- 
evation) in order to increase the exposure time, still keeping a 
reasonably low background. 

We obtained a single spectrum by merging the XIS and 3 
(the two front-illuminated CCDs) data and fitted it simultane- 
ously with the XIS1 spectrum. All the fits were done in the 1-10 
keV range. Both a power law and a blackbody model gave ac- 
ceptable results, with the best fit parameters reported in Table [2] 
These values are similar to those derived in the XMM-Newton 
observation. 

We searched for periodicities in the 2-1000 s range by an 
FFT analysis on the merged event lists from the three instru- 
ments. No periodicities were found, even trying different energy 




274.6 



Fig. 3. Image of the region of AX J1818.8— 1559 in the 2-10 keV en- 
ergy range obtained with the XMM-Newton EPIC instrument. The cir- 
cle (radius 1.8') indicates the error region of the INTEGRAL burst. 
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Fig. 4. Near infrared (K-band) image of the region of AX 
J1818.8— 1559 from the UKIDSS. The circle with radius 0.6 arcsec in- 
dicates the source position derived with Chandra. 



selections. The 3a upper limit on the 1-10 keV pulsed fraction 
is ~30-40% between 2 and 40 s, and ~25-30% above 40 s. 



5. Other archival data 

5.1. Swift 

A short observation (2 ks) with Swift/XRT was performed on 
2007 October 19 about 2 days af ter the INTEGRAL detec- 
tion of the burst dEvans et al.ll2.Q07l) . A source with count rate 
0.012±0.002 counts s _1 was detected at the coordinates R. A. = 
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Table 2. Spectral results. 





XMM-Newton 
2003 


Chandra 
2007 


Suzaku 
2011 


r 

N H (10 22 cm" 2 ) 
F a (10~ 12 erg cm -2 s _1 ) 
F b (10~ 12 erg cm -2 s" 1 ) 
X 2 /dof, dof 


1.17±0.17 
3.6±0.5 
1.37±0.07 

2 11+0.17 

0.88, 79 


1.56 + 0.22 

4 K+0.6 

1 02+ 07 
1.10,91 


1.53 + 0.14 

6 3 +0 ' 8 
1.04 ± 0.03 

2 1Q+ 016 
1.07, 155 


kT BB (keV) 
N H (10 22 cm" 2 ) 


1.87+0.12 
1.6+0.3 


1 42+ - 11 
2.5 + 0.3 


1,<J *— 0.08 
o c+0.5 



R^ B (km) 0.12+0.05 0.16 + 0.02 0.11 + 0.01 

F a (10~ 12 erg cm" 2 s" 1 ) 1.26+0.07 8.4tg j 1.00 ± 0.03 

F b (10 -12 erg cm -2 s" 1 ) 1.45+0.06 1.09 + 0.05 1.21 + 0.04 

X 2 /dof, dof 0.92, 79 1.12,91 1.07,155 



Quoted errors are at la confidence level. 
a 2-10 keV flux. 

b 0.5-10 keV flux corrected for absorption. 
d For d = 10 kpc. 



18 h 18 m 51M9, Dec. = -15°59'20'.'07 (J2000; 4'.'1 error radius, 
90% confidence) consistent with those of AX J1818.8 — 1559. 
The small number of counts is not sufficient to perform a spectral 
analysis. Assuming the parameters of the XMM-Newton power- 
law fit, we derived a 2-10 keV observed flux of —1.5 x 10~ 12 
erg cm~ 2 s _1 . 

This region of sky was observed again with the Swift/XRT 
on 2012 March 04, but no source was detected within or close 
to the error region of the burst. The short exposure time (0.5 ks) 
and off-axis position result in a rather shallow upper limit on the 
count rate of 0.027 counts s _1 (3 a c.l), which is a factor — 2 
above the flux level measured in October 2007. 

5.2. Chandra 

The Chandra observations of the field of AX J1818.8— 1559 
are listed in Table [J0. The 2007 observations were done with 
the ACIS-S instrument as a Target of Opportunity 5 days after 
the INTEGRAL burst detection. They consist of 5 ks of imag- 
ing data, with time resolution of 3.2 s, and a longer exposure in 
continuous clocking (CC) mode, which provides a higher time 
resolution (2.85 ms), but imaging only along a single direction. 
AX J1818.8 — 1559 is detected with a significance above lOer in 
the imaging data, at coordinates R.A. = 18 h 18 m 51 f ?38 (±0?06), 
Dec. = -15°59'22'.'62 (±0'.'05) (J2000); the quoted la uncer- 
tainties are statistical only, and much smaller than the Chandra 
absolute astrometric accuracy, 0'/6 (radial, 90% c.l). We show in 
Fig.|4]the Chandra error circle overlaid on the near infrared (ri- 
band) image from the Galactic Plane Survey, collected as a part 
of the Unit e d Kin gdom Infrared Deep Sky Survey (UKIDSS, 
iLucas et all (120081) ■). No counterparts are visible down to the 
limit K~ 17-18. 

We extracted the spectra from the 2007 observations as fol- 
lows. For the CC data, the source photons were accumulated 
within a 5 x 5 pixels region centred on AX J18 18.8 — 1559 
(one ACIS-S pixel corresponds to 0.492 arcsec); the background 
events were extracted from two symmetric (with respect to the 
source) regions along the image strip. For the imaging data, the 



3 Another observation was performed in August 2009 (obs. ID 
10503), but the position of AX 11818.8-1559 fell in the gap between 
the ACIS-S CCDs 



source counts were selected in a — 5-px radius circular region 
and the background in an annulus with radii of 20 and 30 pixels. 
After checking that they gave consistent results, the two spectra 
were fitted together yielding the best fit parameters reported in 
Table|2] 

AX J1818.8— 1559 was also observed two further times in 
2008. In the first data set (2008 February) the source is de- 
tected at a significance level of 5a, but the small statistics do 
not allow us to carry out a detailed analysis. We just note 
that the 0.3-7 keV count rate, 0.026 ± 0.004 counts s~\ is in 
line with that observed in the 2007 ACIS-S imaging observa- 
tion, 0.024 ± 0.002 counts s _1 . In the second data set (2008 
September), a short HRC-I exposure, the source is not detected. 
The 3a upper limit on the count rate is 0.026 counts s _1 which, 
adopting the Chandra power-law model of Table [2] translates 
into an upper limit on the observed 2-10 keV flux of —5 x 10 -12 
erg cm~ 2 s _1 . 

The most sensitive data for timing analysis are provided by 
the 2007 data in CC mode, which yielded — 1400 source photons. 
By an FFT analysis of these data, we could set a 3a upper limit 
of —40% on the 1-10 keV pulsed fraction for periods in the range 
0.01-1000 s, and —40-60% for periods between 5 and 10 ms. 

5.3. ROSAT 

The position of AX J18 18.8 — 1559 was serendipitously covered 
by a ROSAT observation on 1 1-12 September 1993, falling at an 
off-axis angle of —35' in the PSPC instrument. From these data, 
providing an exposure of —7 ks we derived a 3a upper limit of 
—0.01 counts s _1 on the source count rate in the 0.1-2.4 keV 
energy range. Assuming the same spectral shape of the XMM- 
Newton observation, this corresponds to an unabsorbed flux of 
the order of 2x 10~ n erg cm~ 2 s _1 in the 0.5-10 keV range. 



6. Discussion 

A sub class of GRBs, the so called X-ray Flashes dHeise et all 
12001b . can have spectra as soft as observed in GRB 071017. 
However, in the class of short GRBs with duration less than 2 
s, such soft spectra are never observed. On the other hand, the 
hard X-ray (E>20 keV) spectra of short bursts from SGRs are 
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usually well described by thermal bremsstrahlung models wit h 
temperature -20-40 keV (lAptekar et al .112001 [ iGotz et al.ll2004 . 
as observed in GRB 071017. Assuming a distance of 10 kpc, 
the peak luminosity of the INTEGRAL burst is of the order of 
— 5 x 10 39 erg s -1 , which is a reasonable value for a Galactic 
SGR. Thus, although we cannot completely exclude that GRB 
071017 was a GRB seen through the Galaxy, we believe that, in 
view of its location in the Galactic plane and spectral properties, 
an explanation in terms of a Galactic AXP/SGR is much more 
likely. 

The chance probability of finding an X-ray source of a 
given brightness inside the 10 arcmin 2 INTEGRAL error cir- 
cle can_be_e^tmiated_using the LogN-LogS relations derived 
by ISugizaki et al.l ([2001) from the survey of the Galactic plane 
in which AX J1818.8 — 1559 was discovered. Including the 
small contribution from extragalactic sources seen through the 
Galactic plane, these relations give a surface density of 1.2 
deg~ 2 , for sources with flux equal or greater than that of AX 
J1818.8 — 1559. This leads to a probability of - 3 x to find 
AX Jl 8 18.8 -1559 in the INTEGRAL error region by chance, 
which indicates its very likely association with the burst. 

The X-ray spectrum of AX J1818.8— 1559 is harder than 
that of typical SGRs and AXPs, which, when fitted with a 
single power law, give photon index in the ra nge — 2 — 4 
(Mereghet ti et al.ll2002t iKaspi & Bovdstunl EzOlO). Harder spec- 
tra have been observed in transient magnetars during the initial 
phases of their outbursts, but our analysis of X-ray observations 
of AX Jl 8 18.8 -1559 obtain ed from 2003 to 2012 , as well as the 
earlier ASCA observations (ISugizaki et al.ll200ll) . do not show 
evidence for a transient behavior. Indeed all the measured fluxes, 
and the upper limits, indicate a steady luminosity of 2 x 10 34 
(d/[10 kpc]) 2 erg s~\ consistent with that observed in persistent 
magnetars. The lack of detection of pulsations is not particularly 
constraining, considering that many of the AXPs and SGRs have 
pulsed fractions sm aller than the upp er limits we could set on 
AX J1818.8-1559 dMereghetti et alJ200l . 

Only a single burst has been observed from the direc- 
tion of AX J 1 8 1 8.8 — 1559, despite this region of sky in 
the Galactic bulge having been repeatedly observed by high- 
energy satellites. This contrasts with the abundance of bursts, 
often grouped in time intervals of spasmodic activity, ob- 
served in the most active membe rs of the SGR f amily (e.g., 
SGR 1806-20 (Go giis etal.ll2000h. SGR 1 900+14 (llsrael et all 
120081). SGR 1627-41 dEsposito et al.l l2008l). or IE 1547.0-5408 
dMereghetti et al.l l2009h ). The behavior of AX J1818.8— 1559 
is instead similar to that of a few magnetars from which bursts 
were detected only spora dically, such as for e xample IE 1048.1- 
5937 dGavriil et al.l2002l) . SGR 0418+5729 (Ivan der Horst et all 



I2010I). SGR 1833- 0832 dGogiis et al.ll2010h . or 4U 0142+61 
dGavriil et al.l[20TTh . 

We finally note that the available observations do not favour 
alternative explanations for the nature of AX J1818.8 — 1559. 
Coronal emission from a star can be ruled out by the lack of 
a bright optical counterpart, the hard spectrum, and the large ab- 
sorption, inconsistent with the small distance that would be re- 
quired in this hypothesis. The X-ray to optical flux ratio excludes 
the possibility of a high mass X-ray binary. A low mass X-ray 
binary or an AGN cannot be excluded, but they are not supported 
by the hard spectrum and lack of long term variability, respec- 
tively. 



7. Conclusion 

A new observation of AX J18 18.8- 1559 obtained with the 
Suzaku satellite, as well as a comprehensive analysis of archival 
data from different missions could not establish the real nature 
of this source, which is most likely associated with a SGR-like 
burst observed with INTEGRAL in October 2007. The X-ray 
properties of AX J1818.8— 1559 are consistent with those ob- 
served in AXPs and SGRs, and we propose it as a likely new 
member of this class of sources deserving further study. 
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